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ABSTRACT 

The recent Spitzer detections of the 9.7 fim Si-0 silicate emission in type 1 AGNs pro- 
vide support for the AGN unification scheme. The properties of the silicate dust are of 
key importance to understanding the physical, chemical and evolutionary properties 
of the obscuring dusty torus around AGNs. Compared to that of the Galactic inter- 
stellar medium (ISM), the 10 /im silicate emission profile of type 1 AGNs is broadened 
and has a clear shift of peak position to longer wavelengths. In literature this is gen- 
erally interpreted as an indication of the deviations of the silicate composition, size, 
and degree of crystallization of AGNs from that of the Galactic ISM. In this Letter 
we show that the observed peak shift and profile broadening of the 9.7 /im silicate 
emission feature can be explained in terms of porous composite dust consisting of or- 
dinary interstellar amorphous silicate, amorphous carbon and vacuum. Porous dust is 
naturally expected in the dense circumnuclear region around AGNs, as a consequence 
of grain coagulation. 
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1 INTRODUCTION 

Dust is the cornerstone of the unification theory of active 
galactic nuclei (AGNs). This theory proposes that all AGNs 
are essentially "born equal" - all types of AGNs are sur- 
rounded by an optically thick dust torus and are basically 
the same object but viewed from different lines of sight (see 
e.g. Antonucci 1993; Urry & Padovani 1995): type 1 AGNs, 
which always display broad hydrogen emission lines in the 
optical and have no obvious obscuring effect, are viewed face- 
on which allows a direct view of the central nuclei, while type 

2 AGNs are viewed edge-on with most of the central engine 
and broad line regions being hidden by the obscuring dust. 

Silicate dust, a major solid species in the Galactic inter- 
stellar medium (ISM) as revealed by the strong 9.7 /im and 
18 /im bands (respectively ascribed to the Si-O stretching 
and O-Si-O bending modes in some form of silicate mate- 
rial, e.g. olivine Mg2a;Fe2-2a;Si04), has also been detected in 
AGNs both in emission and in absorption (see Li 2007 for a 
review) . 

The first detection of the silicate absorption feature 
in AGNs was made at 9.7 /im for the prototypical Seyfert 
2 galaxy NGC 1068 (Rieke & Low 1975; Kleinmann et al. 
1976), indicating the presence of a large column of silicate 
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dust in the line-of-sight to the nucleus. It is known now that 
most of the type 2 AGNs display silicate absorption bands 
(e.g. see Roche et al. 1991, Siebenmorgen et al. 2004, Hao et 
al. 2007, Spoon et al. 2007, Roche et al. 2007) which is ex- 
pected from the AGN unified theory - for a centrally heated 
optically thick torus viewed edge-on, the silicate features 
should be in absorption. 

For type 1 AGNs viewed face-on, one would expect to 
see the silicate features in emission since the silicate dust in 
the surface of the inner torus wall will be heated to tempera- 
tures of several hundred kelvin to ~ 1000 K by the radiation 
from the central engine, allowing for a direct detection of 
the 9.7 /im and 18 /im silicate bands emitted from this hot 
dust. However, their detection has only very recently been 
reported in a number of type 1 AGNs covering a broad lu- 
minosity range, thanks to Spitzer (Hao et al. 2005, Sieben- 
morgen et al. 2005, Sturm et al. 2005, Weedman et al. 2005, 
Shi et al. 2006, Schweitzer et al. 2008). It is worth noting 
that the silicate emission features have recently also been 
detected in type 2 QSOs (Sturm et al. 2006, Teplitz et al. 
2006). 

Compared to that of the Galactic ISM, the 9.7 /jm sil- 
icate emission profiles of some AGNs appear "anomalous". 
As illustrated in Figure [1] both quasars (high luminosity 
counterparts of Seyfert 1 galaxies; Hao et al. 2005, Sieben- 
morgen et al. 2005) and the low-luminosity AGN NGC 3998 
(Sturm et al. 2005) exhibit silicate emission peaks at a 
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Figure 1. Comparison of the silicate emission features of the 
quasar 3C 273 (Hao et al. 2005) and the low-luminosity AGN 
NGC 3998 (Sturm et al. 2005) with the silicate absorption feature 
of the ISM toward the Galactic center object Sgr A* (Kemper et 
al. 2004). Most notably is the long wavelength shift of the peak 
positions of 3C 273 and NGC 3998 in comparison with the ISM 
profile. All profiles are normalized to their peak values. 

much longer wavelength (~ 10-11.5 /im), inconsistent with 
the "standard" silicate ISM dust (which peaks at ~9.7/im). 
The 9.7 /im feature of NGC 3998 is also much broader than 
that of the Galactic ISM (Sturm et al. 2005) Q 

The deviations of the silicate emission profiles of type 1 
AGNs from that of the Galactic ISM dust are generally in- 
terpreted as an indication of the differences between AGNs 
and the Galactic ISM in the composition, size distribution, 
and degree of crystallization of the silicate dust (Sturm et 
al. 2005). However, we show in this Letter that the observed 
peak shift and profile broadening of the 9.7 pm silicate emis- 
sion features of AGNs can be explained in terms of porous 
composite dust consisting of ordinary interstellar amorphous 
silicate, amorphous carborQ and vacuum, without invoking 
an exotic dust composition or size distribution. In the dense 
circumnuclear region around AGNs, a porous structure is 
naturally expected for the dust formed through the coagu- 
lation of small silicate and carbonaceous grains. 



2 MODEL OF POROUS COMPOSITE DUST 

To model the porous composite dust, we adopt the multi- 
layered sphere model originally developed by Voshchinnikov 

1 We should note that the width of the interstellar 9.7 /an Si— O 
absorption feature is not universal, but varies from one sightline 
to another (see Draine 2003). Generally speaking, it is relatively 
narrow in diffuse clouds and broad in molecular clouds (Bowey et 
al. 1998). In contrast, its peak wavelength is relatively stable. 

2 There must exist a population of carbonaceous dust in the AGN 
torus, as revealed by the detection of the 3.4 (im absorption fea- 
ture, attributed to the C— H stretching mode in saturated aliphatic 
hydrocarbon dust (see Li 2007 for a review). Whether the bulk 
form of the carbonaceous component in AGNs is amorphous car- 
bon or hydrogenated amorphous carbon is not clear. 



& Mathis (1999). This model assumes that the dust con- 
sists of many concentric spherical layers of different types 
of materials. Each of the material has a given volume 
fraction. Voshchinnikov et al. (2005, 2006) demonstrated 
that the optical properties of porous dust calculated from 
the multi-layered sphere model is in close agreement with 
that from the discrete dipole approximation (DDA; Draine 
1988). While the DDA method is computationally very time- 
consuming, the multi-layered sphere model is computation- 
ally much less demanding but still accurate for computing 
the integral scattering characteristics (e.g. extinction, scat- 
tering, absorption cross sections, albedo and asymmetry pa- 
rameter) and becomes a robust tool to model the optical 
properties of composite dust with a porous structure. 

We consider three types of dust materials in the model: 
amorphous silicate with olivine composition (MgFeSi04), 
amorphous carbon, and vacuum. We take the optical con- 
stants of amorphous olivine MgFeSi04 from Dorschner et al. 
(1995); for amorphous carbon, we take those of Rouleau & 
Martin (1991). 

In the multi-layered sphere model, we need to specify 
(1) flayer - the number of layers; (2) P - the dust poros- 
ity (i.e. the volume fraction of vacuum in a fluffy porous 
grain); (3) m ca rb/m s ii - the mass ratio of amorphous car- 
bon to amorphous silicate in a grain; (4) r compac t - the ra- 
dius of the mass-equivalent compact sphere. Voshchinnikov 
et al. (2005) found that the optical properties of layered 
dust are independent of the number and position of layers 
when the number of layers ni aycr exceeds 15. We therefore 
take ni a y er = 180 We consider a range of porosities, with 
P ranging from P — (compact dust) to P = 0.9. We 
take the mass ratio of amorphous carbon to amorphous sil- 
icate to be m car b/m s ii = 0.7 which is estimated from the 
cosmic abundance constraints (see Li & Lunine 2003)0 We 
take the radius of the mass-equivalent compact spherical 
grain to be r compac t =0.1 pan, a typical size for interstellar 
dust0 For a porous dust of porosity P with the same mass 
as that of the compact dust of radius r compact , its radius is 

^porous — f"comp a ct/(l P) ^ • For a given ^compact, ^porous 

moderately increases with P. 



3 RESULTS 

Using the computational techniques of the multi-layered 
sphere model of Voshchinnikov & Mathis (1999) and as- 
suming niayer = 18, m carb /m si i = 0.7, and r compa ct = 
0.1 pm (see §2), we calculate the absorption efficiency fac- 
tors Qabs(A) for porous composite dust of a range of porosi- 
ties. We show in Figure [2] the 8-13 /jm absorption effi- 
ciency Q a bs(A) calculated for porous composite dust con- 

3 This means that each layered dust has 6 shells, and every shell 
has 3 layers consisting of amorphous silicate, amorphous carbon 
and vacuum. 

4 With the mass density of amorphous silicate p s ;i = 3.5 gem -3 
and the mass density of amorphous carbon p car b = 1.8gcm — 3 , 
this mass ratio corresponds to a volume ratio of V car b/Ksii ~ 1.4. 

5 With a larger r compa ct , the peak wavelength Apg^ and width of 
the 9.7pm Si-O band will be further red-shifted and broadened. 
But as ^compact exceeds ^ 0. 5 /im, the 9.7 /xm Si— O feature fades 
away for porous dust with P > 0.5. 
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Figure 2. 8-13 fim absorption efficiencies Q a bsM of porous 
composite grains consisting of amorphous silicate, amorphous car- 
bon and vacuum with P = 0, 0.2, 0.5, 0.7, 0.9. Note the progressive 
broadening of the 9.7 fim Si-O feature and the progressive shift 
to longer wavelengths of its peak position with the increase of the 
dust porosity P. 

sisting of amorphous silicate, amorphous carbon and vac- 
uum with P = 0,0.2,0.5,0.7,0.9. With r compact = 0.1 ^m, 
the radius of the porous dust corresponds to r polous ~ 
0.108, 0.126, 0.149, 0.215 /im for P = 0.2, 0.5, 0.7, 0.9, respec- 
tively (for compact dust, the 9.7 (im silicate absorption pro- 
files are essentially identical for grains of radii r — 0.1, 0.108, 
0.126, 0.149, 0.215 (jm, as expected since they are in the 
Rayleigh regime at A = 10 /jm). Most notably in Figure[3]are 
the progressive broadening of the 9.7 /jm Si-O feature and 
the progressive shift to longer wavelengths of the peak po- 
sition of this feature as the dust porosity P increases: while 
the silicate feature peaks at ~ 9.6 fim and has a FWHM (full 
width half maximum) ~ 2.1 fim for compact dust (P = 0), 
its peak shifts to ~ 10.6 fim and its FWHM is broadened to 
~ 2.8 /im for very fluffy dust (P = 0.9) ; the 9.7 fim Si-O fea- 
ture is significantly flattened as the porosity increases from 
P = to P = 0.9. 

In order to have a direct comparison between the porous 
dust model with the "anomalous" silicate emission features 
observed in AGNs, we first multiply the calculated absorp- 
tion efficiency Q a bs(A) with a Planck function B\(T) at tem- 
perature T. For a given T we further multiply the product 
of Qabs(A) and B\(T) with a constant to force the model 
to fit the observed flux density of the 9.7 fim feature. This 
approach is valid if the silicate feature emitting regions are 
optically thin0 

6 The optical thin assumption is reasonable if the silicate emis- 
sion mainly comes from the unblocked surface layer of the inner 
torus wall. However, the inferred temperatures (T < 220 K) for 
3C 273 and NGC 3998 (see Figs.Eflit are much lower than that ex- 
pected for the dust near the surface of the inner torus wall where 
the temperature of silicate dust should be close to its sublimation 
temperature (~ 800-1500 K; Kimura et al. 2002). Sturm et al. 
(2005) suggested that this emission may actually originate from 
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Figure 3. Comparison of the silicate emission spectra calcu- 
lated from the porous composite dust model with that observed in 
3C 273 (a bright quasar) and NGC 3998 (a low luminosity galaxy). 
The model emission spectra are obtained by folding the absorp- 
tion efficiencies Q a bsM of porous composite dust (see Fig. [2} with 
a Planck function at temperature T and then scaled to the flux 
densities of the 9.7 fim features of 3C 273 and NGC 3998. 

We take 3C 273 (a bright quasar; Hao et al. 2005) and 
NGC 3998 (a low luminosity AGN whose luminosity is ~4- 
5 orders of magnitude below those of quasars; Sturm et al. 
2005) as two test cases. While their 9.7 fim Si-O features are 
virtually identical (see Fig.[T|, at the 18 fim O-Si-0 feature 
region they deviate significantly from each other (see Fig. [3]). 

As shown in Figure |3 with T~ 120-220 K, the 
porous composite dust model fits the broadened and long 
wavelength-shifted 9.7 ^m feature of 3C 273 and NGC 3998 
reasonably well. The required temperature (T~ 120-220 K) 
is consistent with that constrained by Hao et al. (2005) 
and Sturm et al. (2005). With T = 130 K the model (with 
P < 0.2) is also in a reasonably good agreement with the 
long-wavelength wing of the 18 fim feature of 3C 273. It is 
more challenging to fit that of NGC 3998 which is much 
weaker than that of 3C 273 and peaks at ~ 18.5 fim (while 
the "18 ^m" O-Si-O feature of 3C 273 peaks at ~20/im)Q 

the extended narrow-line regions (NLRs) which are beyond the 
torus. The NLR origin of the silicate emission is consistent with 
the detection of silicate emission in type 2 QSOs (see Efstathiou 

2006, Marshall et al. 2007, Schweitzer et al. 2008). Alternatively, 
the silicate emission could originate from the inner regions of the 
torus provided that the torus is clumpy so that the light from the 
central engine is substantially attenuated (e.g. see Levenson et al. 

2007, Nenkova et al. 2002, 2008). 

7 It is possible that the composition of the silicate dust in 
NGC 3998 may differ from that in 3C 273, suggesting that there 
might be significant environmental variations (after all, NGC 3998 
is a LINER galaxy with an AGN luminosity ~1.7 X 10 4 times 
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Qualitatively speaking, one requires a higher T and a high 
porosity (P > 0.7) to fit the 18 fim feature of NGC 39980 
We should emphasize that the major purpose of this work 
is not to provide a detailed modeling of the silicate emis- 
sion spectra of 3C 273 and NGC 3998, but to put forward a 
hypothesis that the peak redshift and profile broadening of 
the 9.7 /mi Si-O emission features observed in some AGNs 
could result from the porous structure of the dust. 

In the Hao et al. (2005) sample, four quasars (among 
the five PG quasars displaying the 9.7 /im and 18 /im sil- 
icate emission features, apart from 3C 273) exhibit an 
18 fim O-Si-O emission feature intermediate between that 
of NGC 3998 and that of 3C273 (while their 9.7 pm fea- 
tures are all similar). As can be seen in Figure |3j the porous 
composite dust model with T ~ 130-150 K should be able to 
simultaneously fit both the 9.7 /im features and the 18 pm 
features of those quasars. 



4 DISCUSSION 

It is well-known theoretically that as the size of a sili- 
cate grain increases, the 9.7 /mi Si-O feature becomes wider 
with its peak shifted to longer wavelengths (see Fig. 9 of 
Dorschner et al. 1995, Fig. 1 of Voshchinnikov & Henning 
2008). But for compact silicate spheres to have its 9.7 /im 
Si-O feature peaking at A pea k > 10.5 /jm, their size needs to 
exceed ~ 2 /im. However, for these large grains the 9.7 /im 
Si-O feature fades away (see Fig. 6 of Greenberg 1996, Fig. 1 
of Voshchinnikov & Henning 2008). Therefore, it is difficult 
to account for the "anomalous" silicate emission features of 
AGNs just in terms of an increase in grain size. 

A non-spherical grain shape or shape distribution can 
also broaden the 9.7 /jm Si-O feature and red-shift its peak 
wavelength (see Li 2008). But to account for the peak wave- 
lengths of Apeak ~ 10-11.5 /im observed in some AGNs, the 
dust has to be extremely elongated which is unrealistic (e.g. 
spheroidal dust needs to have an elongation >6). 

The peak wavelength of the 9.7 /jm Si-O feature also 
depends on the silicate mineralogy. But we are not aware of 
any amorphous silicate species that could give rise to a Si-O 
feature peaking at A > 10 /im. 

Crystalline olivine has its Si-O feature peaks at ~ 
11.2 jum (see Yamamoto et al. 2008). But this feature is 
too sharp compared to the broad Si-O emission features 
of AGNs. Within the signal-to-noise ratio limits, we see 
no clear evidence for crystalline silicates in 3C 273 and 

below that of the bright quasar 3C 373; it is therefore possible 
that the silicate dust in NGC 3998 subjects to different degrees 
of processing compared to that in 3C273.) As shown experimen- 
tally in Dorschner et al. (1995), the peak wavelength, width and 
strength of the "18 /an" O-Si-O feature (relative to the 9.7 (im 
Si-O feature) vary among silicate minerals of different composi- 
tion. The optical constants of amorphous olivine MgFeSiCU cho- 
sen here may not be the most suitable ones for NGC 3998 (e.g. 
the fast declining red wing of the 18 (im feature may suggest the 
presence of clino-pyroxenes [Wooden et al. 1999]). 
8 One may ask why our models imply that the dust in 3C 273 
is cooler than the dust in NGC 3998 while 3C 273 is much more 
luminous. A plausible answer is that the silicate emission actu- 
ally originates from the NLRs far away from the central heating 
regions or from a clumpy, attenuated torus. 




Porosity 



Figure 4. Peak wavelengths of the "9.7 fim" Si— O silicate fea- 
tures of Q a bsM or QabsM B\(T) as a function of dust porosity 
for a range of temperatures T = 100, 130, 150, 180, 200 and 250 K. 
Also plotted are the peak wavelengths of the 9.7 /im feature of 
NGC 3998 and 3C 273. The shaded area plots the range of the 
peak wavelengths of the 9. 7 /mi Si— O emission features of the five 
quasars observed by Hao et al. (2005): 10 (im < A pca k < 11.5 /im. 



NGC 3998 (see Fig.Q}. So far, the detection of crystalline sili- 
cate dust in AGNs has only been reported in the BAL (broad 
absorption line) quasar PG 2112+059 (Markwick-Kemper et 
al. 2007)13 The 9.7 /im Si-O feature emission features of 
the protoplanetary disks around T Tauri and Herbig Ae/Be 
stars are often also much broader than that of the ISM (e.g. 
sec Bouwrnan et al. 2001, Forrest et al. 2004). This is usually 
interpreted as an indicator of grain growth (e.g. see Natta 
et al. 2007). The porous composite dust model was recently 
used by Voshchinnikov & Henning (2008) to demonstrate 
that a similar behavior of the feature shape occurs when the 
porosity of fluffy dust varies. Kriigel & Siebenmorgen (1994) 
have also demonstrated that a porous structure results in the 
broadening, weakening and redshifting of the 9.7 /im Si-O 
feature. Fluffy dust aggregates are also naturally expected 
in protoplanetary disks as a result of grain coagulation. 

We plot in Figure [4] as a function of porosity P the 
peak wavelengths A pca k of the 9.7 /im Si-O feature of the 
absorption efficiency profiles Q a bs(A) and the emission pro- 
files obtained by folding Q a bs(A) with the Planck function 
B\{T) at various temperatures. It is seen that for 3C 273 and 
NGC 3998, to account for the observed A pca k ~ 10.6 /mi, 
one requires either cool dust (with T ~ 140-150 K) with a 
low porosity (P < 0.3), or warm dust (with T~ 170-220 K) 
with a high porosity (P > 0.6). The 18 /im O-Si-O emis- 



9 Spoon et al. (2006) reported firstly the detection of narrow ab- 
sorption features of crystalline silicates in ultraluminous IR galax- 
ies (ULIRGs). 
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sion feature further constrains that cool dust is preferred in 
3C 273 while NGC3998 appears to favour warm dust. 

Also seen in Figure [4] is that the range of the peak 
wavelengths 10 fim < A pca k < 11.5 )im of the 9.7 /im Si- 
O emission features of the five PG quasars of Hao et al. 
(2005) falls within the predictions of the porous composite 
dust model. This indicates that the long-wavelength shifted 
9.7 /im Si-O features of these quasars can be accounted for 
by the combined effects of porosity and temperature (i.e. 
either highly porous warm dust or cold dust with a lower 
porosity). The 18pm O-Si-O emission feature will allow 
us to break this degeneracy. Finally, for those AGNs whose 
9.7 fj,m Si-O absorption or emission features do not exhibit 
any long- wavelength shift, one may resort to dust with a low 
porosity P < 0.3 (and T > 220 K if they are in emission). 

We should stress that although the combined effects of 
porosity and temperature suffice by themselves to explain 
the general trend of broadening the 9.7 pm Si-O emission 
features of AGNs and shifting their peak wavelengths to 
longer wavelengths, we by no means exclude the effects of 
other factors such as grain size0 shape, and mineralogy. It 
is very likely that all these factors act together to produce 
the anomalous silicate emission features seen in AGNs. 

Finally, we emphasize that although a steeply rising 
(cold) Planck function could redshift the 9.7 /im Si-O emis- 
sion feature, the observed shift of the peak wavelength of this 
emission feature in AGNs cannot be purely a temperature 
effect since the silicate absorption profiles of some AGNs 
also appear anomalous; e.g., the "9. 7pm" silicate feature of 
Mkn 231, a peculiar type 1 Seyfert galaxy, is seen in absorp- 
tion peaking at ~ 10.5 /im (Roche et al. 1983); Jaffe et al. 
(2004) found that the 9.7 fim silicate absorption spectrum of 
NGC 1068 shows a relatively flat profile from 8 to 9 /im and 
then a sharp drop between 9 and 10 /im; in comparison, the 
Galactic silicate absorption profiles begin to drop already at 
~ 8 /mi. 

To conclude, we have explored in a quantitative way 
on the effects of grain porosity on the silicate Si-O stretch- 
ing feature using the multi-layered sphere model. It is found 
that the Si-O feature broadens and shifts to longer wave- 
lengths with the increasing of dust porosity. We conclude 
that the combined effects of dust porosity and cool tem- 
perature of T < 200 K (which further redshifts the silicate 
feature) could explain the observed broadening and longer- 
wavelength shifting of the 9.7 fim Si-O feature of AGNs. 



10 The grain size effect is probably (at least in part) responsi- 
ble for the nondctection of the 9.7 fim and 18 fim silicate emis- 
sion features in many type 1 AGNs (e.g. see Hao et al. 2007). 
The silicate size distribution in some AGNs might be dominated 
by large grains (> a few /im; see Maiolino et al. 2001a, b) or 
small silicate grains are depleted (e.g. see Laor & Draine 1993, 
Granato & Danese 1994). Alternative explanations include so- 
phisticated torus geometries (e.g. tapered disk configurations [Ef- 
stathiou & Rowan- Robinson 1995], clumpy torus models [Ncnkova 
et al. 2002, 2008; but see Dullemond & van Bemmel 2005]) and an 
assumption of strong anisotropy of the source radiation (Manske 
et al. 1998). 
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